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Abstract. We consider the cross section in Fourier space, conjugate to the outgoing hadron's trans- 
verse momentum, where convolutions of transverse momentum dependent parton distribution func- 
tions and fragmentation functions become simple products. Individual asymmetric terms in the cross 
section can be projected out by means of a generalized set of weights involving Bessel functions. 
Advantages of employing these Bessel weights are that they suppress (divergent) contributions from 
high transverse momentum and that soft factors cancel in (Bessel-) weighted asymmetries. Also, 
the resulting compact expressions immediately connect to previous work on evolution equations for 
transverse momentum dependent parton distribution and fragmentation functions and to quantities 
accessible in lattice QCD. Bessel-weighted asymmetries are thus model independent observables 
that augment the description and our understanding of correlations of spin and momentum in nu- 
cleon structure. 
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Introduction 

In the factorized picture of semi-inclusive processes, where the transverse momentum 
of the detected hadron P/,j^ is small compared to the photon virtuality Q^, transverse 
momentum dependent (TMD) parton distribution functions (PDFs) characterize the spin 
and momentum structure of the proton [1, 2, 3, 4, 5, 6, 7]. At leading twist there are 8 
TMD PDFs. They can be studied experimentally by analyzing angular modulations in 
the differential cross section, so called spin and azimuthal asymmetries. These modula- 
tions are a function of the azimuthal angles of the final state hadron momentum about 
the virtual photon direction, as well as that of the target polarization (see Fig. 1 and e.g., 
Ref. [8] for a review). TMD PDFs enter the SIDIS cross section in momentum space 
convoluted with transverse momentum dependent fragmentation functions (TMD FFs). 
However, after a two-dimensional Fourier transform of the cross section with respect to 
the transverse hadron momentum Ph±, these convolutions become simple products of 
functions in Fourier bj-space. The usefulness of Fourier-Bessel transforms in studying 
the factorization as well as the scale dependence of transverse momentum dependent 
cross section has been known for some time [9, 10, 11, 12, 13, 14, 15]. Here we exhibit 
the structure of the cross section in b^-space and demonstrate how this representation 
results in model independent observables which are generalizations of the conventional 



weighted asymmetries [16, 6, 7]. Additionally, we explore the impact that these observ- 
ables have in studying the scale dependence of the SIDIS cross section at small to mod- 
erate transverse momentum where the TMD framework is designed to give a good de- 
scription of the cross section. In particular we study how the so called soft factor cancels 
from these observables. The soft factor [17, 15, 18, 14, 19, 20] is an essential element 
of the cross section that emerges in the proofs of TMD factorization [11, 13, 14, 15]. 
It accounts for the collective effect of soft momentum gluons not associated with either 
the distribution or fragmentation part of the process and it is shown to be universal in 
hard processes [18]. 



Weighted asymmetries 

The concept of transverse momentum weighted Single Spin Asymmetries (SSA) was 
proposed some time ago in Ref. [6, 7]. Using the technique of weighting enables one 
to disentangle in a model independent way the cross sections and asymmetries in terms 
of the transverse momentum moments of TMDs. A comprehensive list of such weights 
was derived in Ref. [7] for semi-inclusive deep inelastic scattering (SIDIS). In SIDIS 
and Drell-Yan scattering, proofs of TMD factorization contain an additional factor, the 
soft factory [17, 15, 18, 14, 19, 20]. At tree level (zeroth order in as) the soft factor is 
unity, which explains its absence in the factorization formalism considered for example 
in Ref. [8]. Consequently it is also absent in tree level phenomenological analyses of 
the experimental data (see for example Refs. [21, 22, 23]). Potentially the results of 
analyses at different energies can be difficult to compare . For a correct description of 
the energy scale dependence of the cross sections and asymmetries involving TMDs, 
the soft factor is essential to include [19, 14, 24]. However, it is possible to consider 
observables where the soft factor is indeed absent or cancels out. These are precisely the 
weighted asymmetries. 

We will focus on the Sivers asymmetry [2]. With a general |P/,x| -weight, this asym- 
metry can be written as [7] 
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The above asymmetry has already been discussed for the choice 

In the numerator, the angular weight sin(^/, — ^5) projects out the structure function 

F^j^P '^^^ [8] from the cross section, while the |P/jx | -weight leads to a "deconvolution" 
of that structure function into a product of the first pr -moment of the Sivers function 

f^j^^ {x) and the lowest moment of the unpolarized fragmentation function D^^^ {z) [7]. 

In [25], we show that in TMD factorization the soft factor cancels in the asymme- 
try, due to the "deconvolution" achieved by appropriate |P;,^|- weighing. We demon- 
strate that it is a natural generalization to employ Bessel functions as weights, w„ oc 
Jn{\^hi-\^T), where 3§t is the Fourier conjugate variable of |P/,_l|. 



Moreover, this generalization addresses a problem related to the perturbative tail of 
TMDs. Without further regularization, the integrals defining the moments f^j^\x) and 
f^^\x) are ill-defined due to the asymptotic behavior /j^(;c,p|-) °^ l/Pj-, fii^^v'r) °^ 
1/Py at large pr (see [26]). Using Bessel functions as weights, the respective integrals 
become convergent while the "deconvolution" property and soft factor cancellation are 
preserved. However, it is important to stress that while the integrals are convergent, the 
scale dependence of the resulting functions, and consequently also the dependence 
of the asymmetries, remains to be studied. Also the link with gluonic pole matrix 
elements [27] becomes less direct. 



Bessel- weighted asymmetries 

To deconvolute or convert the convolutions of TMD PDFs and TMD FFs in the SIDIS 
cross section into products, one can perform a multipole expansion and a subsequent 
Fourier transform of the cross section with respect to the transverse components P/,j^ 
of the hadron momentum. In general, we can write a transverse momentum dependent 
cross section C7(|P/jj^ | , in the form 
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where d'(br) = l^J^=_o„e'"'^''a„(|br|) is a two-dimensional multipole expansion of the 
cross section in Fourier space. The harmonic in 0/, is accompanied by the Bessel 
function of the first kind 7„. 

With these definitions, the relevant terms of the SIDIS cross section can be written as 
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(3) 



where the ellipsis represent 16 more terms. Note that there is only a finite number 
of multipoles in the SIDIS cross section; the Bessel function of highest order is 73. 
Here we show only two terms in the cross section and omit for now regularization 
parameters needed beyond tree level, see Ref. [25] and references therein for more 
details. Introducing the Fourier-transformed TMDs and fragmentation functions 

f{x,hl) ^ J jV^'P^-''^/(^,Pr) = 2;r^°°J|pr||pr|7o(|br||pr|)/(^,P^) , 



FIGURE 1. Kinematics of the SIDIS process. The in- and out-going lepton momenta are / and /', 

respectively. The momentum transfer is q. The target nucleon carries momentum P and its transverse spin 
components are labelled S_l. The momentum of the measured hadron has transverse components Pftj_, 
which define an angle with the lepton plane. 



D(z,b|) = y^/2Kre'''^-^^£>(z,K|)=2;ryd|Kr||Kr|yo(|br||Kr|)D(z,4)(4) 

and the derivative (or br moment) 
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the structure functions in Fourier space in Eq. (3) are given by 

Fuu,T = xl^eUl'^''ix,z'h^T)Df^''iz,hl)S{hl)HuuAQ'), (6) 

a 



where, /j*^^" and f^j^^" are the Fourier transformed unpolarized and Sivers TMD PDFs 

respectively, and D^^^" is a Fourier transformed fragmentation function. We have used 
the kinematic variables = ~q^, = P^, x ^ xb = Q^/P ■ q, y = P ■ q/P ■ I, and 
Zh = P ■ Ph/P ■ Q aiid assume M <^ Q, |P/i_l| -C zQ. The sum Y,a i^i^s over quark 
flavors a and Ca is the corresponding electric charge of the quark. In contrast to the tree- 
level equation [8, 25], we include here an explicit soft factor 5(by) (in Fourier space) and 
a hard part H{Q^), as in reference [19, 14]. For brevity we suppress the dependencies on 
a UV cutoff scale /i and on rapidity cutoff parameters (e.g., ^, t,, p in [19, 14, 25]). It 
is now clear that the cross section equation (3) is a multipole series, and that projection 
on Fourier modes in polar coordinates (j)h, \^h±\ will give access to the right hand sides 
of equations (6) and (7). Calculating the weighted asymmetry equation (1) with weights 
wi{\Fh±\) = 2Ji{\Fh±\^T)/zM^T and wo(|P/,±|) = M\Ph±\^T) thus yields 

A^^ [X,y,Z,^T) - ^j:„,2^,,^,(g2)/(0)«(,,,2^2),-(^2)5(0)«(^^^,) •(«) 

Due to the "deconvolution" of the structure functions in the weighted asymmetries, and 
universality of the soft factor, S cancels in the numerator and the denominator. 

Further, note that 3§t enters the weights wq and as a free parameter that we can 
scan over a whole range in order to compare the transverse momentum dependence of 



the distributions in the numerator and denominator relative to each other (in Fourier 
space). At the operator level, ^t{— |br|) controls the space-like distance between 
quark fields in the correlation functions we measure. The Bessel-weighted asymmetries 
are a natural extension of conventional weighted asymmetries [6, 7] with weights wi 
proportional to powers of |P;,^ | . Indeed, in the limit 0, equation (8) results in the 

often encountered special case for the SIDIS cross section 



(9) 



In particular, we obtain formally the relations (x, 0) = (x), f^j^^ (x, 0) = f^j^^ ix) 
and (x, 0) = (x), where 

/W(x,0) =/W(x) ^ I d^T (||2)"/(x,p|). (10) 

An analogous equation holds for the fragmentation functions D. However, we caution 
the reader that these moments are not well-defined without some regularization. It is 
therefore safer to study Bessel-weighted asymmetries at finite SSj, where the Bessel 
functions suppress contributions from large transverse momenta. 



Cancellation of soft factor at the level of matrix elements 



In a similar manner to the discussion above, we now consider the soft factor cancel- 
lation in the average transverse momentum shift of unpolarized quarks in a transversely 
polarized nucleon for a given longitudinal momentum fraction x. This shift is considered 
in [28] and defined by a ratio of the pr -weighted correlator. 



S±=0,Sr=(l,0) 



(11) 



where f^j^^ and are the moments defined in Eqs. (10). Obviously, the average mo- 
mentum shift is very similar in structure to the weighted asymmetry Eq. (8). While 
the weighted asymmetries are accessible directly from the P;,j^- weighted cross section, 
the average transverse momentum shifts are obtained from the pr -weighted correlator 
and could in principle be accessible from weighted jet asymmetries. As already men- 
tioned, the integrals defining the moments of TMD PDFs on the right hand side of the 
above equation are divergent without suitable regularization. Therefore we generalize 
the above quantity, weighting with Bessel functions of Iprl. In particular, we replace 



Py = |pr| sin(</)p) 



2/i(|pr|^r) 
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where <^s = for the choice Sj- = (1,0) in Eq. (11). The correlator reads in 

terms of the amplitudes a|^^ and B^'^^ [25], 
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The Bessel-weighted analog of Eq. (11) is thus 
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Again, the soft factors cancel where the independence of the soft factor on v-b/vv"^ is 
crucial [25]. Further, weighting with Bessel functions at various lengths thus allows 
us to map out, e.g., ratios of Fourier-transformed TMDs [25]. In the limit 0, 
we recover equation (11), {pyix))ju = {py{x))Tu > which we have thus shown to be 
formally free of any soft factor contribution. However, we caution the reader again that 
the expressions at — can be ill-defined without an additional regularization step. 



Conclusions 

We have demonstrated that rewriting the SIDIS cross-section in coordinate space dis- 
plays the important feature that structure functions become simple products of Fourier 
transformed TMD PDFs and FFs, or derivatives thereof. The angular structure of the 
cross section naturally suggests weighting with Bessel functions in order to project 
out these Fourier-Bessel transformed distributions, which serve as well-defined replace- 
ments of the transverse moments entering conventional weighted asymmetries. In addi- 
tion, Bessel-weighted asymmetries provide a unique opportunity to study nucleon struc- 
ture in a model independent way due to the absence of the soft factor which as we have 
shown cancels from these observables. This cancellation is based on the fact that the soft 
factor is flavor blind in hard processes, and it depends only on by , /i^, p. Moreover, evo- 
lution equations for the distributions are typically calculated in terms of the (derivatives 
of) Fourier transformed TMD PDFs and FFs. As a result the study of the scale de- 
pendence of Bessel-weighted asymmetries should prove more straightforward. For the 
above stated reasons we propose Bessel-weighted asymmetries as clean observables to 
study the scale dependence of TMD PDFs and FFs at existing (HERMES, COMPASS, 
JLab) and future facilities (Electron Ion Collider, JLab 12 GeV). 
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